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Low-lying one-particle anomalous excitations are studied for Gutzwiller-projected strongly cor- 
related BCS states. It is found that the one-particle anomalous excitations are highly coherent, 
and the numerically calculated spectrum can be reproduced quantitatively by a renormalized BCS 
theory, thus strongly indicating that the nature of low-lying excitations described by the projected 
BCS states is essentially understood within a renormalized Bogoliubov quasi-particle picture. This 
finding resembles the well-known fact that a Gutzwiller-projected Fermi gas is a Fermi liquid. The 
present results are consistent with numerically exact calculations of the two-dimensional t-J model 
as well as recent photoemission experiments on high-Tc cuprate superconductors. 
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Since their discovery in 1986 [lj], high-Tc cuprate su- 
perconductors and related strongly correlated electronic 
systems have become one of the largest and most impor- 
tant fields in current condensed matter physics Q. In 
spite of enormous theoretical and experimental efforts de- 
voted since their discovery, understanding the mechanism 
of the superconductivity and the unusual normal state 
properties of the high-Tc cuprate superconductors re- 
mains unresolved and there is still extensive debate. This 
is mainly because of the strong-correlation nature of the 
problems, which is widely believed to be a key ingredient. 
Immediately after the discovery, Anderson [3| proposed 
a Gutzwiller-projected BCS state to incorporate strong- 
correlation effects in the superconducting state. While 
the ground state properties of projected BCS states have 
been studied extensively Q , it is only very recently that 
the low-energy excitations of projected BCS states have 
been explored by numerically exact variational Monte 
Carlo techniques [f| H, 0, Q, which in fact have found 
many qualitative as well as very often quantitative sim- 
ilarities to the main features observed experimentally in 
high-Tc cuprate superconductors. It is therefore very im- 
portant and timely to understand the nature of low-lying 
excitations described by Gutzwiller-projected BCS states 
because of their great relevance to the high-Tc cuprate 
superconductors. Furthermore, a Gutzwiller projected 
single-particle state is one of the most widely used corre- 
lated many- body wave functions in a variety of research 
fields [H [l(| [llf, and thus better understanding of the 
nature of projected BCS states is also highly desirable. 
This is precisely the main purpose of this study. 

In this paper, one-particle anomalous excitations are 
studied for Gutzwiller-projected BCS superconductors. 
It is shown that one of the characteristic properties for 
the Gutzwiller-projected BCS states is their highly coher- 
ent one-particle anomalous excitations. It is found that 
the numerically calculated one-particle anomalous excita- 



tions can be reproduced quantitatively by a renormalized 
BCS theory. This finding thus strongly indicates that 
the low-lying excitations of the Gutzwiller-projected BCS 
states are described within a renormalized Bogoliubov 
quasi-particle picture, which resembles the well known 
fact that a Gutzwiller-projected Fermi gas is a Fermi liq- 
uid [Hj]. The present results also provide a theoretical 
justification for utilizing a simple mean-field based BCS 
theory to analyze low-energy experimental observations 
for the high-Tc cuprate superconductors. 

A Gutzwiller-projected BCS state 1*^) witn N elec " 
trons is defined by 
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where |BCS) = riker7ko-|0) is the ground state of the 
BCS mean-field Hamiltonian [l3| with singlet pairing and 
7kcr is the standard Bogoliubov quasi-particle annihila- 
tion operator with momentum k and spin er(=f, |), 
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Vg — ~~ r " l iT"-u) i s the Gutzwiller projection oper- 

ator excluding sites doubly occupied by electrons, and 
Vn the projection operator onto even number N of elec- 
trons. Ck<j = 5Dj e~ lk ' x c.\ a j\fL (L: number of sites) is 
the Fourier transform of the electron annihilation oper- 
ator Cio- at site i with spin a, and fi\ a — -Ci a . Note 

that, since the number N of electrons is even, IvPq^) is a 
spin singlet with zero total momentum. In the following, 
it is implicitly assumed that the gap function in |BCS) 
is real and the spatial dimensionality is two dimensional 
(2D). However, the generalization of the present study is 
straightforward. 

A single-hole (single-electron) added excited state 
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by 
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which has momentum k, total spin 1/2, and z-component 
of total spin a. Hereafter, the normalized wave functions 
for the N- and (JV± l)-particle states are denoted simply 
by \ipo) and l4o- )> respectively. The quasi-particle 
weights for the one-particle added and removed normal 
excitations are thus defined as 



and 
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respectively, where a is the opposite spin to a. 

The one-particle anomalous excitation spectrum is 
generally defined as 



F(k,w) = 



1. 
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where \&n) is the nth eigenstate in = 0, 1, 2, • • • , with 
corresponding to the ground state) of a system described 

by Hamiltonian H with its eigenvalue and N elec- 
trons [TtJ ■ The spectral representation thus leads 
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i.e., as will be shown later, the quasi-particle weight for 
the one-particle anomalous excitations is (k). It is 
also interesting to notice that the above equation relates 
Z [ 2 N \k) to the quasi-particle weights for the one-particle 
normal excitations, i.e., 
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(k) = Zf^-k j)-Zi JV ^(kT), (9) 
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which should be useful for computing Z_ (kcr) [3]. 
The validity of Eq. © can be checked numerically by 



computing all quantities, Z+(— k |), Zi JV+ZJ (k f), and 



/(^+2), 



r(N) 



(k) 



independently. A typical set of results calcu- 
lated by a standard variational Monte Carlo technique on 
finite clusters is shown in Figs[T] (a) and (b), where one 
can see that indeed Eq. © holds within the statistical 
error 12011. 
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Note that F(k, w) is a real function provided that spin 
rotation and reflection invariance as well as time reversal 
symmetry are assumed The frequency integral of 



the one-particle anomalous excitation spectrum, 



(N) 
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(JV) ), 



provides a 



well-known sum rule which will be used later. 

To study the one-particle anomalous excitations for the 
projected BCS state, let us first defined the following 
quantity similar to -F^ for the projected BCS states: 



rf)( k) = 



-kT C -ki 
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FIG. 1: (Color online) (a) and (b): ^ J (k) (cir- 
cles) and Z ( / } (-k |) • Z ( _ N+2) (k T) (crosses) for differ- 
ent momenta k and N = 222 (n w 0.87). (c): I*] 2 = 
KV'o + l c iT c i+x|IV'o JV ')| 2 (circles) as a function of electron 
density n. \ipQ ) is optimized for the 2D t-t'-J model with 
i'/t — —0.2 and J/t = 0.3 on an L = 16 x 16 cluster with peri- 
odic boundary conditions 0, [2lj]. For comparison, zfjI^BCsl 2 
for a renormalized BCS theory (see the text) is also presented 
by crosses in (c). 
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is constructed in exactly the same way as 



\tpQ y> ) except that the number of electrons onto which 
the state is projected is N + 2. Using the previously 
derived relations for the projected BCS states one 
can now easily show that 
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Now let us assume that there exists a system (with 
Hamiltonian H) for which the ground state and the low- 
energy excited states can be described approximately 
by the projected BCS states \^ N) ) and l*^ 13 ), i.e., 
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mediate consequence of this assumption is that the sum 
rule for F(k,w) is now duF{k, w) = Z { 2 N \k). An- 
other consequence is that the spectral representation of 
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the one-particle anomalous excitations is 



F(k, W ) = Zf(k).i( u 
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i.e., the quasi-particle weight for the one-particle anoma- 



lous excitations is Z. 



Here E { Q N) = 

= £(<r } ). 



>(k). 

(4 N) \H\4 N) ), and E { _1 +1) = E(V^>). This is be- 
cause of the equality derived here in Eq. ([5]). Since the 
spectral weight is not positive definite, the above equa- 
tion along with the sum rule does not immediately im- 
ply that the contribution of incoherent "other terms" 
in Eq. (|10|l is negligible. However, using the prop- 
erty of the projected BCS states reported previously 
that the one-particle added normal excitations are co- 
herent ID, 1(1 [22| , one can easily show that indeed the 



one-particle anomalous excitations consist of a single co- 
herent part for each k with no incoherent contributions. 
It is interesting to note that numerically exact diagonal- 
ization studies of small clusters have also found that the 
one-particle anomalous excitations for the 2D t-J model 
are highly coherent with relatively small incoherent con- 
tributions 
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Let us now calculate the one-particle anomalous excita- 
tion spectrum F(k, uj) for the projected BCS states. For 
this purpose, here we will consider the 2D t-t'-J model 
on the square lattice @. This model has been studied 
extensively and found to show a d-wave superconducting 
regime in the phase diagram [It], Hj|- Furthermore, it is 
well-known that a Gutzwiller-projected BCS state with 
rf-wave pairing symmetry [Eq. |T])] is a faithful variational 
ansatz for the superconducting state of this model 
The model parameters used here are set to be t' jt = —0.2 
and J/t = 0.3 [24 1 . 

The results of the numerically calculated F(k, u>) for 
representative momenta are shown in Fig. [21 where 
is optimized to minimize the variational energy 
= 222 on an L = 16 x 16 cluster (n w 0.87) with 



2l|. As is expected for 



for N 

periodic boundary conditions 
a c?-wave superconductor, the spectral weight becomes 
smaller toward the nodal line in the (0, 0)-(7r, 7t) direc- 
tion [see also Figs. Q] (a) and (b)], and it changes the sign 
across the nodal line where the weight is zero. 

To understand the nature of the low-lying excitations 
observed in F(k,w), here the results are analyzed based 
on a renormalized BCS theory with d-wsve pairing sym- 
metry [13j. The procedure adopted is as follows: (i) the 

E { n N) is fitted for 



(N+l) 



excitation energy E(k) = E, 
all momenta k in the whole Brillouin zone by a standard 
Bogoliubov excitation spectrum [25| , (ii) using the fitting 
parameters determined in (i), the BCS spectral weight for 

the one-particle anomalous excitations, m^ 80 ^ t>£. BCS L is 
calculated, and (iii) the BCS spectrum is renormalized 
by a momentum independent constant zb in such a way 
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FIG. 2: (Color online) The one-particle anomalous excita- 
tion spectrum F(k,uj) for the projected BCS states [Eq. (I10|l ] 

(solid lines). ) use d here is optimized for the 2D t-t'- 

J model with t' /t = -0.2, J/t = 0.3, and N = 222 on an 
L = 16 x 16 cluster (n ~ 0.87) with periodic boundary con- 
ditions 21]. The momenta k studied are indicated in the 
figures. For comparison, F(k, u>) for a renormalized BCS the- 
ory (see the text) is also presented by dashed lines. The 
momentum-independent renormalization factor zb is 0.30. 
The delta function 5(ui) is represented by a Lorentzian func- 
tion e/n(Lu 2 + e 2 ) with e = 0.02i. 
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The obtained renormalized BCS spectra are shown in 
Fig.[5]by dashed lines. It is clearly seen in Fig.[5]that the 
renormalized BCS spectra can reproduce almost quanti- 
tatively F(k, uj) for the projected BCS states. It should 
be emphasized that the procedure employed above is 
highly nontrivial and it is beyond a simple fitting of nu- 
merical data. Similar agreement is also found for differ- 
ent sets of model parameters, one of which is exemplified 
in Fig. [3] The surprisingly excellent agreement found 
here strongly indicates that the low-lying excitations de- 
scribed by the projected BCS states can be well under- 
stood within a renormalized Bogoliubov quasi-particle 
picture. 

To further examine the validity of the renormal- 
ized Bogoliubov quasi-particle picture for the projected 
BCS states, let us finally study the superconduct- 
ing order parameter, which is here defined as \& = 

(V'o Ar+2 L c i| c i+xil'^o Ar ' ) ) ( x being the unit vector in the 
x direction). The electron density (n) dependence of W 
is shown in Fig. [T] (c) for the 2D t-t'-J model, where 
\ipQ ) is optimized for each n 21|. As seen in Fig. [1] 
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FIG. 3: (Color online) The same as in Fig.[2]but for N = 240 
(n ~ 0.94) [2lJ. The momentum-independent renormalization 
factor zb is 0.18. 



(c), 1^1 2 vs n shows a domelike behavior, similar to the 
pairing correlation function at the maximum distance 
as a function of n reported before ||. It is also inter- 
esting to notice that |\&| 2 is proportional to 1 — n for 
small 1 — n. The corresponding quantity ^bcs for the 
BCS state with uj^ ^ and ^ BCS \ determined by the 
procedure mentioned above, can also be calculated by 



*bcs = iEk e ""I w k ■ If a renormalized Bo- 
goliubov quasi-particle picture is valid, 'F = zs^bcs is 
expected. As seen in Fig. [T] (c), this is in fact clearly the 
case. This result also gives a clear physical meaning to 
the renormalization factor zb introduced in Eq. pip . 

As is well known, a Gutzwiller-projected Fermi gas 
is described within a Fermi liquid picture The 
present results thus strongly suggest that analogously a 
Gutzwiller-projected, correlated BCS state ("projected 
BCS gas") can still be described within a renormal- 
ized BCS-Bogoliubov quasi-particle picture ("BCS liq- 
uid" ) [2^| . This is in fact in accordance with recent pho- 
toemission spectroscopy experiments on high-Tc cuprate 
superconductors for which low-lying excitations consis- 
tent with a BCS theory have been revealed j22|- More- 
over, the present results would also provide a theoreti- 
cal justification for employing a mean-field-based BCS- 
like theory to analyze the low-energy dynamics observed 
experimentally in the superconducting state of high-Tc 
cuprate superconductors [28}. 

To summarize, the one-particle anomalous excitations 
have been studied to understand the nature of the low- 
lying excitations of strongly correlated superconductors 
described by the Gutzwiller-projected BCS states. It was 
found that the low- lying excitations, which are highly 



coherent, can be essentially described within a renor- 
malized Bogoliubov quasi-particle picture. This finding 
thus resembles the well-known result that a Gutzwiller- 
projected Fermi gas is a Fermi liquid. Finally, the present 
study has demonstrated that a variational Monte Carlo- 
based approach can be also utilized to explore low-lying 
excitations, and hopefully this work will stimulate further 
studies in this direction for other dynamical quantities. 
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